Trapping Reactions with Randomly Moving Traps: Exact Asymptotic Results for 

Compact Exploration. 



G.Oshanin^, O.Bcnichou^, M.Coppey^, and M.Moreau"'^ 
^ Laboratoire de Physique Theorique des Liquides (CNRS - UMR 7600), Universite Pierre et Marie Curie, 4 place Jussieu, 

75252 Paris Cedex 05, France 

^ Laboratoire de Physique de la Matiere Condensee, College de France, 11 place M.Berthelot, 75231 Paris Cedex 05, France 

(Received: June 17, 2002) 

In a recent Letter Bray and Blythe have shown that the survival probability PAit) of an A particle 
diffusing with a diffusion coefficient Da in a ID system with diffusive traps B is independent of Da 
in the asymptotic limit t —* oo and coincides with the survival probability of an immobile target 
in the presence of diffusive traps. Here we show that this remarkable behavior has a more general 
range of validity and holds for systems of an arbitrary dimension d, integer or fractal, provided that 
the traps are "compactly exploring" the space, i.e. the "fractal" dimension d^ of traps' trajectories 
is greater than d. For the marginal case when d^ = d, as exemplified here by conventional diffusion 
in 2D systems, the decay form is determined up to a numerical factor in the characteristic decay 
time. 
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Trapping A+B ^ B and recombination A+B re- 
actions (TR and RR) involving randomly moving A and 
B particles which react "when they meet" at a certain 
distance b are ubiquitous in nature. A few stray examples 
include quenching of delocalized excitations, coagulation, 
recombination of radicals, charge carriers or defects, or 
biological processes related to population survival . 

In recent years there has been much interest in the 
long-time behavior of these processes, following a remark- 
able discovery |J|,|,|,|[0,|,|,|ig of many-particle effects, 
which induce essential departures from the convention- 
ally expected behavior 

A pronounced deviation from the text-book predic- 
tions was found for the diffusion-controlled RR in case 
when initially the particles of the A and B species are 
all distributed at random with strictly equal mean den- 
sities no- It has been first shown Q and subsequently 
proven that as i ^ oo the mean density n{t) follows 

n{t) ~ adfiy^ {Dt)^'^/'^ , where d is the space dimension- 
ality, ad is a constant and D = Da -f I^b is the sum of 
particles' diffusion coefficients. This law contradicts the 
decay law obtained within the Smoluchowski approach 
(SA): n{t) ~ ll4>i^\t) §, where, as t 



oo. 



ht\t)^ I drKsir) 



^yjDt/TT, d = 1, 

ln(4Dt/fe2) ' = ^' 
A^nDbt, d = 3, 
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Ks{t) being the d-dimensional Smoluchowski- type "con- 
stant" , defined as the flux of diffusive particles through 
the surface of an immobile sphere of radius 6. 

For the TR two situations were most thoroughly stud- 
ied: the case when As diffuse while Bs are static, 
and the situation in which the As are immobile while 
Bs diffuse - the so-called target annihilation problem 
(TAP). In the case of static, randomly placed (with 



mean density p) traps the A particle survival probability 
PAit) shows a non-trivial, fluctuation-induced behavior 
g||i,|,|,|0] 
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which is intimately related to many fundamenal problems 
of statistical physics . 

Survival probability Ptarget{t) of an immobile target 
A of radius h in presence of point-like diffusive traps B 
(TAP) can be calculated exactly for any d (see Refs. jlj] 
and 



Ptarget{t) = CXp - p(t>l\tfj 



(3) 



where (t>lf\t) obeys Eq.(|l|) with Da — 0. Decay forms 



in systems with hard-core interactions between Bs 
or with fluctuating chemical activity [ p^ have been also 
derived. 

On contrary, the physically most important case of TR 
when both As and Bs diffuse was not solved exactly. It 
has been proven ji) that here PA{t) obeys 



\nPA{t) = -Xd{DA,DB) X 



t 

Wr 

t, d 



d= 1, 

d = 2, 
= 3, 



(4) 



which equation defines its time-dependence exactly. On 
the other hand, the factor Xd{DA, Db) remained as yet 
an unknown function of the particles' diffusivities and 
d. Since the time-dependence of the function on the rhs 
of Eq. (Q) follows precisely the behavior of /* drKs (t) , 
one might expect that the SA provides quite an accurate 
description for this situation and following its spirit to set 
Da — supposing that traps diffuse with the diffusion 
coefficient D = Db+Da- As a matter of fact, it has been 
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often tacitly assumed that when both of species diffuse 
PA{t) obeys Eq.(|) with ^['^^(t) defined by Eq.(0) with 
D — Da+Db- On the other hand, it has been shown that 
Xd{DA, Db) is less than the corresponding prefactor in 
Ks{t) and that it may be bounded by a non-analytic 
function of Da and Db jl^l- A perturbative approach for 
calculation of Xd{DA, Db), as well as corrections to the 
SA in ID systems were presented [Q. It has been also 
noticed that Xd{DA, Db) is not a function of D ^ Da + 
Db only, since the diffusion-reaction equation arc not 
separable . This lack of knowledge of the precise form 
of Xd{DA, Db), of course, constitutes an annoying gap in 
the general understanding of the fluctuation phenomena 
in chemical kinetics. 

In a recent Letter Bray and Blythe (BB) made a 
considerable step towards the solution of this general case 
by showing that, surprisingly, the survival probability of 
an A particle diffusing in a ID system with diffusive traps 
is independent of Da in the asymptotic limit t oo and 
coincides with the survival probability of an immobile 
target in the presence of diffusive traps, Eq.(^). The 
convergence to this asymptotic result might be, however, 
rather slow as shows the comparison |lq] with extensive 
numerical simulations |p^ . 

One may, however, pose quite a legitimate question 
whether such a remarkable result is constrained to physi- 
cally not very realistic ID systems with conventional dif- 
fusion or if it is just a particular case of a more general 
behavior which persists also for higher dl 

In this Letter we show that indeed this remarkable re- 
sult holds for a more general case. Namely, we show that 
for systems of an arbitrary dimension d, integer or frac- 
tal, the large-t behavior of the survival probability of a 
randomly moving A particle in the presence of randomly 
moving traps is given exactly by the solution of the ex- 
actly solvable TAP, provided that the traps are "com- 
pactly exploring" the space; in other words, the "frac- 

/ B) 

tal" dimension dw of traps' trajectories is greater than 
d For lattice random walks, this corresponds to sit- 
uations in which random walks are recurrent pl| . Ran- 

( B) — 

dom motion with rfu, > 2 is widespread in nature and 
is most often encountered in porous and disordered sys- 
tems, amorphous and polymer materials for which 
systems it will take place in two and even three dimen- 
sions (see, e.g., Refs. [^IH and ||2^). Finally, we exam- 
ine the behavior in the marginal case when dw — d, as 
exemplified here by conventional diffusion in 2D systems, 
and show that here the decay form can be determined up 
to a numerical factor in the characteristic decay time. 

Consider a d-dimensional volume V containing a single 
mobile A particle of radius b and N point-like traps B. 
Let be the vector denoting the A particle position at 
(7) 

time moment t, while , j = 1, . . . ,N , be the corre- 
sponding vector denoting the position of the J-th trap. 
Introducing two auxiliary indicator functions 



<5b(x) 



1, |x|<6, 
0, otherwise; 



i{v) 



1, y = o, 

0, otherwise. 



(5) 



one writes PA{t) down formally as follows: 



„■ 1 



)dT 



}, (6) 



where the symbol E{. . .} denotes averaging with respect 
to the A particle trajectories, while the brackets with the 
subscript {xj"'"'} stand for averaging with respect to the 
trajectories of the the j-th trap. Note that Eq.(||) applies 
for any type of motion provided that the point-like traps 
are ignorant of each other and thus move independently. 
In the limit N,V ^ 00 {N/V = p) one has 

PA{t)^E\e^p(-p f dxo(/'(Xt,Xt)) )|, 

(7) 

where brackets denote now averaging with respect to the 
trajectories of a single trap B whose starting point is at 
position xq, while /'(X(,X() is the indicator function 



l'{Xt,^t) = l-l[J 6b{Xr-Xr)dT 



(8) 



which shows whether two given realizations of trajecto- 
ries Xt and Xf have "intersected" each other at least 
once within the time interval [0,t]. Note that averag- 
ing over trajectories Xi in the exponential is to be taken 
for fixed Xf, and after performing such an averaging we 
have to average an exponential of the result with respect 
to the trajectories Xf, which represents a fairly complex 
mathematical problem. Such a complexity emphasises, 
of course, the significance of the BB result. 

Now, BB have noticed although not proven rigor- 
ously, that the A particle will on average survive longer if 
it stays still than if it diffuses. They have also furnished 
some arguments in favor of this statement showing that 
this is true for systems with a finite number of traps since 
here the lowest value of the decay exponent corresponds 
to Da = 0. In other words, it means that PA{t) in Eq.(0) 
is bounded from above by 



PAit) < Ptargetit) = CXp ( - p4'^\t) 



(9) 



where (t) is given by 



'it) 
= / dxo 



dxo 



I- I 



S],{yLr)dT 



{xt},xt=o=xo 

(10) 



{xt},Xf=o=0 



Note also that Eq. (|9[) should hold for any type of random 
motion, not necessarily only for conventional diffusion. 
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A lower bound on PA{t) can be constructed in the fol- 
lowing way [Q: One notices first that all terms in the 
product in Eq.(||) are positive definite and hence, if one 
performs averaging of Eq.(^) not over all possible real- 
izations of trajectories X( and x^^'-' , but only over some 
restricted subset, one arrives at a lower bound on 
Following Ref. we define this subset as follows: let 
us suppose that initially the A particle has been located 
at the origin while all traps were uniformly spread in a 
d-dimensional system such that the trap nearest to the 
origin appeared at distance I from it. Then, we per- 
form averaging only over such trajectories of the A par- 
ticles which do not leave within the time interval [0, t] 
the volume of radius I centered around the origin, and 
such trajectories of traps B (which all initially are uni- 
formly distributed outside of this volume) do not enter 
there within the time interval [0, t]. For such trajectories 

Y[l[ / <5fc(X.-xW)dr) ^1, (11) 

and hence, the following lower bound is valid 

PAit) > expi-Vdpl"^) X Prob(maa;{|X^|} < 1\t G [0,t]) 
X Probj(mm{|x(.^)|} > 1\t G [0,t]), (12) 

where Vd denotes the volume of a d-dimensional sphere 
of a unit radius, while two other multipliers stand for the 
probability that the A particle does not leave a sphere of 
radius I within the time interval [0, t] and the probability 
that neither of traps, initially uniformly distributed with 
mean density p outside this sphere, enters this sphere up 
to time t. Note that exactly the same lower bound has 
been already proposed in Refs. HJ] and 

Now, let us suppose that the mean-square displace- 
ment of the A particle obeys < X^ {DAt)'^^'^'-^\ 
while the MSD of traps follows < x^ >~ {D etf 
d\f^ and diP^ being the "fractal" dimensions of the A 
particle and traps trajectories, respectively. For conven- 
tional diffusion of has d[f^ — d[^^ = 2 for any d. Under 
quite general conditions, for DAt ^ Z'''^ ' the probability 
Prob(TOaa;{|XT-|} < 1\t E [0,t]) can be estimated as pl| 

Prob(maa;{|X^|} < 1\t £ [0,<]) ^ exp{-Pd{DAt)/l'^^^^^ ), 

(13) 

where (3d is a constant dependent on the type of random 
motion and d. On the other hand, one readily notices 
that Probj(mirt{|xr"'^|} > 1\t G [0,t]) is just the proba- 
bility that an immobile target of radius I survives until 
time t in the presence of randomly moving traps, i.e. 

Probj (mm{ I x*.-'^ I } > 1\t £ [0,t]) = exp - p</>|''\t)) , 

(14) 



where 4>i {t) is defined by Eqs.(p^ with h replaced by /. 

We turn next to the most delicate point of our analysis. 
We note first that the definition in the first line in Eq. ( [To| ) 
allows to express (j}f'\t), in virtue of the Gauss theorem, 
as a time-integral of Ks(t) (see Ref. |^). On the other 
hand, the definition in the second line in Eq.(|l^) shows 
that (l)['^\t) can be thought of as the mean volume swept 
by randomly moving fictitious particle of radius I during 
time t, i.e. the mean volume of the so-called "Wiener 
sausage" (see, e.g. Ref. ||l^). Its lattice counterpart is 
known as the mean number of distinct sites visited pl| . 
General properties of such a volume for different types 
of random motion have been first analysed in the pio- 
neering papers by de Gennes [^0|| , in which he studied 
RR involving polymerized particles and TR on the per- 
colation cluster. As well, behavior of 4>f\t) have been 
discussed at length in Ref. within the context of the 
polymer-free voids distribution in polymer solutions. It 
has been shown that depending on the relation between 
dui and d, two completely different types of behavior 
may be observed. The first type of behavior occurs when 
dL^^ < d. In this case (j)f\t), called by de Gennes as 
the "exploration volume", is smaller than the volume x^ 
where the particle is confined. This case is called the case 
of non-compact exploration and here (jif^ {t) ^ ^dt, where 
the prefactor ^d is some function dependent on d and the 
type of random motion. In this case jd is proportional to 
some positive power of /! For lattice random walks this 
regime takes place when random walks are non-recurrent 
PH . In the opposite case when dL^"* > d the behav- 
ior is completely different. Here, the exploration volume 
4>f\t) increases sublinearly with time, the trajectories 
are spatially more confined and most of space inside the 
volume x^ is indeed visited. This case is called the case of 
"compact exploration" (recurrent random walks for the 
lattice counterparts) and here 

<jyf\t) ^ {\^t\f ^ {DBtY'^'"\ t^oo (15) 

What is most important in this case is that the prefactor 
in this asymptotic law is independent of l\ Note also that 
such a behavior is compatible with the Alexander-Orbach 
result 4>f~'^^\t) ~ t'^/Z'^if' for anomalous random walk 
on fractal lattices of dimension d/ p^ ]. 

Now, we note that the function on the rhs of Eq.(|T^) 
is valid for any value of / and consequently, the "best" 
lower bound would correspond to such I which provides 
its maximal value. Focussing next solely on the case of 
random motion with compact exploration, we note that 
in this case the leading large-t behavior of (t>'l'^\t) is in- 
dependent of /, and hence, we have to maximize only the 
product of the first two terms. This yields 

PA{t) > exp ( - a'a p^-^DAty) exp ( - pcf>f\t)) 

(16) 
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where a'^ is a constant, z = d/{d + d[f^) and the asymp- 



totic behavior of 4>f\t) is defined in Eq.(f5) 



For compact exploration we have that d < di, . On 
the other hand, on comparing the growth rate in the ex- 
ponent in the first term on the rhs of Eq. (^6|) against the 
growth rate of (j)f\t) defined in Eq.(^5|), we infer that 
the second multipher determines the overall decay in case 



when d, d'lf^ and d[^^ 



obey: d < dlj < d + d, 



which 



reduces to simple condition of compact exploration for 



d\^\ Further on, for d, d^^ and diP'' which obey 



the double side inequality, we have evidently that the 
leading terms in Eqs.(^ and ( p^ coincide, since 0^''^ (t) is 
asymptotically independent of b and (/)|''-' (t) does not de- 
pend on I. We infer thus that in this quite general case ex- 
act asymptotic solution for trapping reactions A+B — > B 
in which both species move randomly is given by the so- 
lution of the corresponding immobile target annihilation 
problem, which represents a substantial generalization of 
the BB result 

Finally, we analyze the behavior in the marginal case 



using as an example conventional diffu- 
sion in 2D systems. In this case, an asymptotical behav- 



of. 



a(2) 



{t) is well-known [2^ 
offln(4i?B</P)-27 



1 



(ln(4DBt/?2)-27 



[3] 

[4] 

[5] 

[6] 

[7] 
[8] 

[9; 

■[10 



(17)[ll 



where Ai w 0.423 and 7 is the Eulcr constant. On the 
other hand, for standard diffusive motion one has that 
Prob(maa;{|X^|} < 1\t £ [0,t]) obeys Eq.(|3t) in which 

one sets dL"^' = 2 and (32 is the square of the first zero 
of the Bessel function Jq{x). Now, since (j)] (t) is only 
weakly (logarithmically) dependent on I, one may assume 
that the value of I which maximizes the lower bound for 
DAt ^ P is still determine d by the derivative of the first 
two terms on the rhs in Eq.(12), i.e. I (/32-DaV^2/o)^^'*- 
Such an estimate shows then that in 2D systems with 
conventional diffusion the A particle survival probability 
is bounded by: 



1 



Ai 



\n{ADBt/b^) 



+ C'(ln-2(t)) < 



< \u[llPA{t) 
'ln(4//32) 



ln(4£>Bt/6^) 
ATrDstp 

h In(l/2P6') 



< 2- 



\n{DB/DA) 



ln{4DBt/b^) 



+ 0(ln-2(t)) 



[12 

[13 

[14 
[15 

[16 
[17: 
[18 

[19 
[20; 
21 
22 
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Hence, in this marginal case the suitably extended up- 
per and lower bounds determine the decay form up to a 
numerical factor in the characteristic decay time. 

The authors acknowledge helpful discussions with Pro- ^24 
fessors Bray and Blythe on the matters of this paper. 
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